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Abstract. The ACDM model is the most commonly admitted to describe our Universe. In spite of a great 
success with regard to the large scale structure formation, some problems are still unresolved at galactic 
\ scales. Alternative scenarios have to be explored such as modified gravity. We have developed an N-body 

code able to solve in a self consistent way the galactic dynamics in MOND. The first version of the code 
I consists in solving the modified Poisson equation on a uniform Cartesian grid to derive the gravitational force 

O ^. on each particle. With it, we study the evolution of isolated galaxies, like the bar instability, the angular 

I ■ momentum transfer, etc. Galaxies in MOND are found to form stronger bars, faster than in Newtonian 

2 I dynamics with dark matter. In a second step, we implement an adaptive mesh refinement technique in 

the code, allowing to run more contrasted simulations on larger scales, like interacting galaxies. During 
an interaction, the dynamical friction forces are less important in MOND, and merging times are longer 
than in DM models. The different morphologies of interacting galaxies in the two models are discussed. All 
simulations are performed in both frameworks of modified gravity and Newtonian gravity with dark matter 
with equivalent initial conditions. 



> 

, 1 Introduction 

m 

CO ' MOND is an alternative theory of gravitation proposed by Milgrom in 1983, to account for galactic dynam- 
. ics without dark matter. It is based on a modification of the Newtonian gravity below a universal critical 
' acceleration ao ~ 2 x 10~^^ m s~^. In MOND, the Poisson equation is transformed into: 

r- 

O ■ V[/i(|V^|/ao)V$] = 4^Gp, (1.1) 

has the property to be equal to unity for large accelerations (Newtonian regime) and tends to x for the 
. low accelerations (MOND regime). Because of the non-linearity of the MOND gravity, new potential solvers 
^ \ have to be developed for the modified Poisson equation (Brada & Milgrom 1999, Ciotti et al. 2006, Tiret & 
Combes 2007, hereafter TC07). Numerical simulations of spherical dissipationless systems have been carried 
out by Nipoti et al. 2007. They find that the merging time-scales are significantly longer in MOND compared 
to Newtonian gravity with dark matter. 

We present here N-body simulations of isolated galaxies and interacting galaxies with MOND. Model galaxies 
have a stellar and gaseous component. 

2 Numerical technique 

2.1 Uniform grid 

We use a uniform Cartesian grid in three dimensions to study the evolution of an isolated galaxy. The modified 
Poisson equation is solved using full multigrid technique (Tiret & Combes 2007). It is discretized as follows : 

4^Gp,,,, = (2.1) 

(^i+l,j,fc — ^i,j,/c)MMi — (^i,j,/c — ^i-l,j,fc)MLi 
+ (0i,j+l,/c — ^i,j,/c)MM2 — (^i,j,/c ~ ^i,j-l,/c)ML2 
+ ((/>i,j,fe+l - 4^iJ,k)l^M3 — {4^iJ,k — (/>i,j,/c-l)ML3)/^^ 
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Fig. 1. Full multigrid (FMG) algorithm is used to accelerate the convergence in the resolution of the modified Poisson 
equation. 

with Pi,j,/c and (/>i,j,/c the spatial density and potential discretized on a grid of step /i, jUMn f^L^ value of 
jii{x) at points Mi and Li (Fig. [T]). The gradient component {d/dx^d/dy^d/dz)^ in are approximated by 

(/)(/)+0(J/)-(/)(K)-0(J) 0(C)+0(D)-0(£;)-0(F) x 

Our models evolve in a simulation box of 50 kpc^, including 256^ nodes. The spatial resolution is 400 pc 
everywhere. 

2.2 Adaptive refinement resolution 

To simulate interacting galaxies, we now construct an adaptive grid from the mass distribution ([2|) and solve 
the modified Poisson equation (Tiret & Combes in prep). This allows to use a larger simulation box (150 kpc) 
keeping a good resolution (~ 500 kpc) with reasonable CPU time. 

3 Isolated galaxies 

Galaxies are initially built with a Miyamoto-Nagai disc for the star component, a Toomre disc for the gas and a 
Plummer sphere for the bulge. Each galaxy is simulated in the MOND (MOND model) and Newtonian gravity 
with dark matter (DM model). In the DM model, the DM halo is chosen to obtain the same total rotation 
curve than in MOND. 

3.1 Bar formation, resonance 

In TC07, we analyze pure stellar discs. Galaxies in MOND are found to form stronger bars and faster than 
in the DM model. Gas dissipation has now been added using the sticky particle scheme. If the fraction of gas 
is enough (~ 7%, in late type galaxies), it makes the galaxy more unstable to bars. As expected, in the new 
simulations including the gas component, galaxies form bars sooner than without gas, especially in the DM 
model where galaxies are stabilized by the DM halo. 

In Newtonian dynamic with DM, the bar is slowed down by dynamical friction against DM particles. On the 
contrary, in MOND, the bar pattern speed remains constant. It implies differences on the resonance locations 
(like the corotation, ILR, OLR,...). If the bar pattern speed is constant (MOND), the resonance occurs always 
at the same radius. If the bar slows down (DM), the resonances are shifted outward in the disc. Fig. [3] shows 
the apparition of resonant rings and pseudo rings in MOND, corresponding well to the observations. 
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Fig. 2. Adaptive mesh structure during a galaxy merger. 
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Fig. 3. Several examples showing the morphological structures of NGC 2665, ESQ 509-98, UGC 12646 and NGC 1543 
(top panel) compared to simulated galaxies in MOND (bottom panel). Rings and pseudo-rings structures are well 
reproduced with modified gravity. 

3.2 External field effect: EFE 

In a deep MOND regime, the potential of an isolated distribution of finite mass is logarithmic. However the 
internal potential return to a Keplerian behaviour when an external field is applied. It is why an escape velocity 
exists also in MOND if the system is embedded in an external potential (Wu et al. 2007). 

Let us consider a uniform acceleration field. In Newtonian dynamic, a galaxy will be accelerated uniformly 
but the evolution of the galaxy in its rest frame is the same as if there was no external field. Due to the 
non-linearity of the MOND gravity, a galaxy embedded in a uniform acceleration field will have an internal 
dynamics which is different from the isolated case. Gravitational torques act on the galaxy and make it precess 
(Fig. H. 
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Fig. 4. Edge on view (t = 0, 0.5, 1 Gyr ) of a galaxy embedded in an acceleration field tilts from 45° from the galactic 
plane. 



■1., ^ 



Fig. 5. Simulations of the Antennae galaxies in the DM model (left) and MOND model(right). 
4 Interaction of galaxies 

With the adaptive multigrid code, we are now able to study interacting galaxies. 

4.1 The antennae 

Our first aim was to reproduce the antennae galaxies, by comparing the morphology of this system in the 
different models, the result is plotted on fig. [5l 

MOND can form quite extended tidal tails as in the DM model. However, the conditions to obtain a merger 
are noticeably different. In DM, the halo plays an important role to decrease the distance between the two 
galaxies by dynamical friction. If the galaxies start on circular orbits, they will decay and merge in ~ 3 relative 
orbits. In MOND, the dynamical friction is efficient only when the two baryonic galaxies overlap each other. 
Starting from circular orbits, galaxies in MOND could stay bound along a Hubble time without merging. On 
the other hand, if the orbits are elliptical (with an impact parameter lower than the visible size of the galaxy) 
merging is possible in a few orbits too. 

4.2 Tidal dwarf 

Several observations of interacting galaxies show the formation of tidal dwarves (10^, 10^ ^o) at the extremity 
of the tidal arms (ARP105, NGC 7253, ...). In numerical simulations, it is easy to form clumps all along the 
tidal tails, but this dwarves are difficult to obtain with standard initial conditions. Bournaud et al. (2003) 
have shown that if the DM halo is extended at least to 10 times the optical radius, the gravitational potential 
is strong enough to form dwarf galaxies at the tip of the tail. The logarithmic potential of the asymptotic 
behaviour of MOND plays naturally the role of an extended DM halo. Therefore these tidal dwarves are also 
obtained in MOND as can be seen on fig. [6l 
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Fig. 6. Tidal dwarf formation at the tip of the tidal tail in MOND. 



5 Conclusion 

Numerical simulations with MOND become more and more realistic by taking into account gas dissipation and 
star formation. With the adaptive mesh refinement code, we are now able to study isolated galaxies as well as 
interacting galaxies in very short CPU time (several hours of mono-processor). 

One particular effect observed in all these simulations is the lack of dynamical friction experienced by the 
baryonic particles against the DM halo, which does not exist in MOND. For isolated galaxies, the dynamical 
friction of the DM model makes the bar slow down. One way to increase again the bar pattern speed is to 
accrete external gas. For interacting galaxies, dynamical friction against the DM halo makes the galaxies merge 
in less than a Gyr. On the contrary, to obtain a merger in MOND, galaxies require a smaller impact parameter. 

Then large differences are expected for the cosmic merger history in MOND and DM model. A given 
merger will have a longer life-time, even in the appearance of the most perturbed morphology (tidal tails, 
asymmetries...). To discriminate between the two models, it will be necessary to determine how many mergers 
a galaxy experienced in its life, still an uncertain parameter. 

The fact that MOND implies longer merger times can solve the problem of the existence of compact group 
of galaxies. The standard CDM scenario predicts that compact groups should have merged a long time ago 
while too many of them are observed now. 
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